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It is believed that fermions in adjoint representation on single site lattice will restore the center 
symmetry, which is a crucial requirement for the volume independence of large-N lattice gauge 
theories. We present a perturbative analysis which supports the assumption for overlap fermions, 
but shows that center symmetry is broken for naive fermions. 
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1. Introduction 

Eguchi and Kawai proposed [1] that a pure gauge S\J(N) lattice theory can be reduced to a 
single site at the limit N — > °°. The proof was based on the assumption that the single site has 
a center symmetry associated with the traces of the Polyakov loops. However, later the center 
symmetry was shown to be broken in the continuum limit [2]. 

It was immediately clear that EK-reduction might work for theories with supersymmetry [3]. 
This has been explored recently in the continuum [4], where it was suggested that fermions obeying 
periodic boundary conditions in the adjoint representation, QCD(Adj), would be volume indepen- 
dent. A continuum perturbation theory analysis on R 3 x S l showed Z# symmetry to be restored. 
Similar calculations on S 3 x S 1 also showed a center symmetric phase [5]. A perturbation theory 
analysis on lattice with Wilson fermions and with one small direction agreed with the result [6, 7] 
(see also [8]). The perturbation theory results have also been confirmed by numerical studies on a 
single site lattice [9]. 

In this work we describe a perturbative calculation of QCD(Adj) on a single site lattice with 
naive and overlap fermions. The results show that the center symmetry is not restored for naive 
fermions but it is restored for overlap fermions for a wide range of Wilson mass values m. The 
study of naive fermions is motivated since we do not expect doublers on a single site lattice with 
adjoint fermions. The complete study including lattice simulations is presented in [10]. 



2. Eguchi-Kawai model with adjoint fermions 

Eguchi-Kawai model with fermions is defined with the action 

S = S g + S f , (2.1) 

where 

4 

S* = -bN £ Tr [UpUvUlut] (2.2) 

and 

S f = -/logdet#„ j0 = -/Trlog# n , 0) (2.3) 

where Up £ SU(/V) are the link matrices, b = ^ is the gauge coupling constant, H n o is the naive 
or overlap fermion operator, and / is the number of Dirac fermion flavors. 

Matrices are the link matrices in adjoint representation calculated from as 

V ab = -Tr[T a UT b U\ (2.4) 

where T a are the hermitian generators of SU(N) normalized as Tr[T a T b ] = 2S ab . 

Both H n and H a are 4(N 2 — 1) x 4(N 2 — 1) hermitian matrices and correspond to naive Dirac 
fermions and overlap Dirac fermions respectively. The exact forms are given in [10]. The determi- 
nant of H n>0 is positive definite and therefore the logarithm is well defined. Furthermore, we can 
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define a hermitian matrix £ 

( n \ 

(2.5) 

(2.6) 
such that, 

I.H n>0 Z = H* i0 (2.7) 

which implies that all eigenvalues of H n o are doubly degenerate reflecting the adjoint nature of 
the fermions. In addition, both naive and overlap fermions obey chiral symmetry and therefore the 
eigenvalues of H n o will come in ± pairs. Therefore the factor, /, in front of can be an integer 
(single Dirac flavor) or half-integer (single Majorana flavor) for all values of N. 1 

3. Weak coupling expansion 

Our aim is to find out if the symmetries are broken in the weak coupling limit. We follow [2] 
and perform the weak coupling analysis by decomposing according to 

U ll = i a *D lt e- u *\ D^e^Sij. (3.1) 

Keeping 6'^ fixed, we expand in powers of a^. The lowest contribution to S g comes from the 
quadratic term in [2] and the lowest contribution to 5/ comes from setting = 0. 

The computation of the fermion determinant reduces to a free field calculation at this order 
and the result is 

S n ,o = -4/ £ In A,, o (0'' - + ) - 4(jV - 1 )/ln A,, (</> ) , (3.2) 

¥J 

where e 1 ^ is the boundary condition in the \i direction. The eigenvalues, ±X(p), are two fold 
degenerate and the explicit expressions are given in [10]. The complete result from fermions and 
gauge fields is 




4/lnA„, o (0'-0^ + 0) \-4{N- I) f\nX, h0 {<$>). (3.3) 



Independent of the actual values of 6'^ , the fermion eigenvalues will have (N — 1 ) zero modes 
with periodic boundary conditions when the mass is set to zero. If all the 6'^ are different for each 

jU, then the fermions should not have exact zero modes when we set equal to (o^ — OfA with 
i ^ j. If the fermion spectrum has more than (N — I) zero modes, we will refer to the extra ones as 
doubler zero modes. 

In order to find the minimum of S, we consider the Hamiltonian 

H = ^{4) 2 + PS. (3.4) 



'Note that one should have written YsH no in (2.3) but this is the same as writing H no as long as / is an integer 
multiple of i . 
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Figure 1: Naive fermions: Left: as a function of / at N = 23, j3 = 4, and ji = 0.01. Right: Smallest 
eigenvalue, X{p) with Pfl = (8^-9^ and i ^ ./', for N = 23, j3 = 4, and = 0.01. 



For large j3, the Boltzmann measure e~ H will be dominated by the minimum of S. We can perform 
a HMC update of the n, 6 system to find this minimum. 

A choice for the order parameters associated with the symmetries is [2] 



1 ( 1 



^Isin 2 !^-^ 



(3.5) 



The value = i corresponds to unbroken symmetry in ;U direction and P^ = totally broken 
symmetry. 

3.1 Naive fermions break the Zjy symmetries 

We assume periodic boundary conditions. We pick one value of N and /3 and calculate P^ 
as a function of /. In order to clearly see symmetry breaking we use rotational symmetry on the 
lattice and choose to label our directions such that Pi < Pi < P3 < P4 for each configuration in 
our thermalized ensemble. We set pL = 0.01 to avoid potential singularities that could occur for the 
massless case. The plot for j8 = 4 confirms the breaking of the center symmetry left panel of Fig. 1 . 
The value is well below \ . 

With N = 23, we expect 22 exact zero modes for A(p) as explained in section 3. We plot the 
average of the smallest eigenvalue, A (p) with = (^6'^ — 6 J ^j and i 7^ j, in right panel of Fig. 1. It 
is clear that this eigenvalue is non-zero for all values of / indicating that there are no doubler zero 
modes. The reason for the breaking of the Zjy symmetries can be understood by looking at the total 
action obtained from (3.3): 



I in 

¥J 



I>in : 



1 



2/£ln 

¥J 



M 2 +I> 2 (e<-e£ 



(3.6) 



The fermionic contribution cannot separate 6'^ = from 6'^ = 6^ + 71 implying that the fermion 
contribution alone will result in a distribution of eigenvalues restricted to a width of %. 
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Figure 2: Overlap fermions: Left: Plot of P^ as a function of / at N = 23, j3 = 1, and n = 0.01. Right:Plot 
of Pfi as a function of malN = 23, /3 = 1, and fi = 0.01 for two different values of /. 



3.2 Overlap fermions do not break the Zjy symmetries 

In contrast to naive fermions the overlap fermions do not break the Zjy symmetries. A plot of 
for several values of / at N = 23 and j3 = 1 with /I = 0.01 in right panel of Fig. 2 shows this to 
be the case. The very small deviation close to / = I is a consequence of finite iV effects. Since we 
probably are not able to go beyond N = 23 in the full simulation of the model, this plot will serve 
as a guide to what one can expect in a full simulation. 

Because there are no doubler zero modes, we do not have any restriction on the values for 
the Wilson mass, m, used in the Wilson-Dirac kernel. However, we cannot make it arbitrarily 
large since one can see by a direct computation that the large m limit of overlap fermions is naive 
fermions [11]. A plot of P^ as a function of m is shown form / = A and / = 1 in right panel of 
Fig. 2. It indicates that 3 < m < 8 is an appropriate range of values of m where the symmetries 
are not broken for / = A and that range only gets bigger as / increases. A plot of the lowest positive 
eigenvalue of H w in left panel of Fig. 3 and the smallest eigenvalue, X (p) with = (^9^ — 6^, 
i 7^ j, in right panel of Fig. 3 shows that there are no doubler zero modes in this range of m. This 
range of m can be used for full numerical simulation with overlap fermions. 

3.3 Effect of fermion boundary conditions 

We have assumed periodic boundary conditions for fermions in the previous two subsections. 
Other choices of boundary conditions that do not generate a U(l) anomaly amount to replacing 
by Vfie 1 f with integer valued [12]. Physical results are expected to depend on the choice of 
boundary conditions. This is in contrast to the case of large ;V gauge theories coupled to fundamen- 
tal fermions. In that case, the e'~N~ factor can be absorbed by a change of gauge fields that only 
changes the Polyakov loop and not the action. If the symmetries are not broken as is the case in 
the confined phase, this change will not affect physical results. 
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Figure 3: Overlap fermions: Left Plot of the lowest positive eigenvalue of H as a function of msiN = 23, 
j8 = 1, and jX = 0.01 for two different values of /. Right:Plot of smallest eigenvalue, X(p) with pn — 

{% - and i ^ for N = 23, J3 = 1, and fl = 0.01. 



In order to study the effect of boundary conditions on symmetry breaking, we focus on 
f = \ and / = 1. We set 0^ = for n =2,3,4 and varied 0i by setting it to equal to ^= with k 
an integer in the range < k < N/2. The plot of as a function of ty\ is shown in Fig. 4. We 
see that the Zpj symmetry in the pL = 1 direction is broken if fy\ > f . This seems to be the case in 
the limit of large ,/V and seems to be roughly independent of /. Furthermore, the Zn symmetries in 
the other three directions with periodic boundary conditions are not broken. This result could help 
us force feed momentum for quarks in the adjoint representation. In order to pursue this, we need 
to study the effect of </>i on the chiral condensate and see if chiral symmetry is restored when the 
symmetry is broken and if the chiral condensate is independent of the value of 0i when the Z^ 
symmetry is not broken. 



4. Discussion 

We have provided perturbation theory arguments that the Zjy symmetry is restored in the large 
N limit for overlap fermions but not for naive fermions. The details of perturbation theory calcu- 
lations with framework to HMC-lattice simulations are given in [10]. The natural continuation is 
to perform lattice calculations on several N and calculate physical observables, e.g., string tension 
and chiral condensate. 
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Figure 4: Overlap fermions: Plot of P^ as a function of <j>\ for three different values of (f,N) with j3 = 1 
and n= 0.01. 
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